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Abstract 

Structure of pipe networks minimizing the total energy dissipation rate is studied analytically. Among all the pipe networks 
possible to build with a given amount of total pipe volume (or pipe lateral surface area) , the network minimizing the dissipation 
rate is shown to be loopless. Furthermore, such an optimal network is shown to contain at most N — 2 nodes, in addition to 
the N sources plus sinks that it connects. These results are valid whether the possible locations for the additional nodes are 
chosen freely or from a set of nodes (such as points of a grid). Applications of these results to various physical situations and 
to efficient computation of optimal pipe networks are also discussed. 

PACS numbers: 89.75.Fb, 05.65.+b, 45.70.Vn, 45.70.Qj 



Finding the most efficient transport network is an issue 
arising in a wide array of different contexts Q, Q ■ One 
can cite, among others, the water, gas and power supply 
of a city, telecommunication networks, railway or car traf- 
fic, and more recently the design of labs-on-chips or mi- 
crofluidic devices. Besides, this problem also appears in 
theoretical works intending to describe the architecture 
of vascular systems of organisms 0, Q . Generally speak- 
ing, consider a set of sources and sinks embedded in a 
two- or three-dimensional space, their respective number 
and locations being fixed. Flow rates coming in through 
each source and going out through each sink are also 
given. The challenge consists in interconnecting them via 
possible intermediate junctions, referred to as additional 
nodes, in the most efficient way. That is, to minimize a 
cost function of general form ^ fc Wkf (ik), where the sum 
is carried on all the links that compose the network. Wk is 
the "weight" associated with the fcth link, and / is some 
function of the flow rate ik carried by this link (usually, 
/ is some power law: / (ik) — |«a| 7 , with 7 > 0). The 
flow rates are not independent but must satisfy a conser- 
vation law (Kirchhoff's first law) at every source, sink, 
and additional node of the network. 

In the present letter, the structure of pipe networks 
minimizing the dissipation rate U = ~^2 k r^i? is studied, 
where the weight rk of each pipe is the " flow resistance" , 
defined as: 
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p being some positive constant, Ik and Sk the respective 
length and cross-sectional area of each pipe, and m a pa- 
rameter characterizing the flow profile (examples of flows 
are given hereafter) . It is shown that - among all the pipe 
networks possible to build with a given value of total pipe 
volume (or total lateral surface area) - the network min- 
imizing U is necessarily loopless. As a consequence, it is 
also shown that the number of additional nodes in such 
an optimal network cannot exceed N — 2, where N is 



the number of initial nodes (sources plus sinks). These 
results are valid whether the positions of the additional 
nodes are chosen in a set or free to be adjusted for a 
further minimization of U. 

Since the minimization is done under a global con- 
straint on the total pipe volume, or pipe surface area, 
it is convenient to introduce the function C n = J2k ^ kS k > 
which represents the total volume (n = 1) or total lateral 
surface area (n = 1/2) of the network jj]. For most flows 
encountered in physics, m > 1, ensuring that m > n. 
In many physical situations also, ik derives from a po- 
tential function (electrical potential, pressure, concentra- 
tion, temperature,...) so that potential difference Vk, flow 
rate ik, and resistance of the fcth pipe are related by 
the Ohm's law Vk — rkik- 

Consider a network of given topology, connecting the 
sources to the sinks, and whose total pipe volume/surface 
area is equal to C n . The dissipation rate U is a function 
of the cross-sections and lengths of all the pipes that com- 
pose the network. Pipe cross-sections can be adjusted to 
minimize the dissipation rate, while preservingthe value 
of C n . It has been shown in a previous study |j| that, as 
a result of this optimization of the network geometry, the 
flow rate ik carried by each pipe in the optimal network 
scales with its cross-sectional area Sk as: 
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where I is the total flow rate coming in through the 
sources (and going out through the sinks), and k is a 
function depending on m, n, and the geometry and topol- 
ogy of the network. Consequently, the dissipation rate in 
such a network can be rewritten as: 



U = (kI) C n 



(3) 



Similarly, pipe lengths may also be optimized in order 
to minimize U (while preserving C n ). Actually, coordi- 
nates of the additional nodes are the appropriate inde- 
pendent optimization parameters. As a result, the fol- 
lowing vector balance is also satisfied at every junction 
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of the network with optimized cross-sections and node 
locations Ul: 
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where is the outward-pointing unit vector along each 
adjoining pipe. In some situations however, the positions 
of the additional nodes cannot be freely adjusted, but 
must be chosen from a set of nodes (such as points of 
a grid, or some particular cities of a country). In these 
cases Eq. 01 is usually not satisfied. 

Conditions and 0] do not uniquely determine the 
structure of the optimal network: they do not restrict 
the number of pipes and junctions, or give indication on 
the presence or absence of loops in the network. Besides, 
these are properly necessary conditions for extrema of U. 
It can be shown (see Appendix) that these extrema are 
either local minima, or saddle points. Here, it is proven 
that condition [3 cannot lead to a local minimum of U 
in a network containing loops. As a consequence, it is 
shown that, starting from any network containing loops, 
it is always possible to build a new loopless network with 
lower dissipation rate (and with a same value of C n ). 

Consider a network satisfying Eq. and containing a 
loop. To go from a given junction A to another junction 
B of the loop, there are two different paths, noted (a) 
and (/?) , as depicted on Fig. QJi. Let us make a shift of 
material, in such a way that flows in path (a) tend to be 
strengthened in one direction (say A to B) and flows in 
path (/?) tend to strengthened in the opposite direction 
(B to A). That is, the new cross-sectional areas s' k in 

the loop are defined as: s 
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with, for path (a), a plus sign if flow rate in pipe is 
in direction A — > B and a minus sign if the flow rate is in 
opposite direction, while signs are inverted for path (/3) 
(see Fig. $0 is a positive number smaller than any 

cross-sectional area Sk of the original loop. Cross-sections 
outside the loop remain unaltered (s' k = Note that 
flows in a loop cannot turn all clockwise, or counter- 
clockwise (otherwise, the potential difference Va — Vg 
between nodes A and B, and thus the flow rates in the 
loop, would trivially be zero). This guarantees that the 
cross-sectional areas of the loop did not all simultane- 
ously increase (or decrease). 

Such a variation of cross-sectional areas implies a re- 
distribution of flows in the entire network. Let {i k } be 
the new distribution of flow rates in this new geometry, 
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the new resistances, and U' = 
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new dissipation rate. Although we do not know the val- 
ues of the new flow rates, an upper bound on the new 
dissipation rate U' can be established, using Thomson's 
principle 0, ■ This principle (which happens to be a 
theorem) states that - among all possible flow rate dis- 
tributions {jk} which satisfy conservation of flow at every 
source, sink, and additional nodes - the actual flow rate 




(a) 



FIG. 1: Shift of material in a loop of the network, (a): the 
original loop, where flow directions in each pipe are indicated 
with arrows, (b): the same loop, where cross-sectional area 
of a pipe is increased when the direction of its carried flow 
is A — > B along path (a) or B —* A along path (/?), and 
decreased otherwise. The other cross-sectional areas in the 
network remain unaltered. 



distribution (i.e.: the one satisfying Ohm's law) is the 
one that makes the function ^2 k Tkj\ an absolute mini- 
mum. Consider in particular the flow rate distribution 
defined as: jk = ik + io along path (a), jk = ik — io along 
path (/3), and jk = ik for any pipe outside the loop. i$ 
is some positive number, and {ik} is the actual distribu- 
tion in the original network, the sign of ik being defined 
in both paths as positive if directed from A to B. The 
distribution {jk} satisfies flow rate conservation at every 
source, sink, and additional nodes, since the distribution 
{ik} does. Besides, by choosing i = KlsQ n+n '' 2 and us- 
ing Eq. [21 the flow rate distribution jk can be rewritten: 
Jk = sgn (i k )Kls k {m+n)/2 . Thus, according to Thomson's 
principle: 



U' <{KlfC n , 



(•5) 



with C" n = y^ fe Ik s'k ■ Let us now compare the new value 
of pipe volume/surface area C' n with the original value 
C n . This can be done by studying the variation of C' n 
with so- With no difficulty, we obtain the derivative of 
this function with respect to x = s^ 71+n ^ 2 : 
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the Since m > n, is a decreasing function of x, and then: 
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Therefore, if the bound in inequality [7| is negative, C' n 
is a decreasing function of x. Using Eqs ^ an d [3 this 
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bound can be rewritten as: 

v 7 ^=0 \path a path ] 

(8) 

We could have chosen to reinforce flows in direction 
B — > A in path (a), and A — > £> in path (/?) instead, 
which comes to swapping (a) and in the calculations 
above. Inequalities|51and|3would still be satisfied for this 
new geometry, but this time with an opposite sign for 

(^ L ) ( see necessarn y (^&) — ^ f° r 

one of the two shifts, and C' n is a decreasing function of s 
for this particular shift, implying that C' n [sq) < C' n (0) = 
C n . From Eqs. |3]and[5l we obtain that the corresponding 
dissipation rate U' is lower too: U' (sq) < U ||. 

Afterwards, the total volume/surface area can be in- 
creased up to its original value C n by increasing any 
cross-sectional areas in the network. This will imply a 
further decrease in U 0. Thus, we found a small pertur- 
bation of the cross-sections such that the dissipation is 
reduced for a fixed value of C n . This proves that condi- 
tion [21 cannot lead to a local minimum of U in a network 
containing loops. Moreover, we can apply the reasoning 
above with increasingly large values of sq, until eventu- 
ally one of the pipe in the loop has a zero cross-sectional 
area, so one of the paths is cut of. Possible dead branches 
can be removed, the equivalent material being shifted to 
the rest of the network by increasing any other cross- 
sectional areas again, so that the constraint stays at its 
initial value while the dissipation rate is subjected to a 
further decrease. We can repeat the same procedure and 
eliminate all the duplicate paths until there are no loops 
in the network. The argument holds even in case of over- 
lapping loops (that is, loops having pipes in common), 
and more generally for any topology of the original net- 
work. Therefore, we conclude that the architecture of 
the network minimizing U is necessarily loopless. Note 
that we do not make use of condition 01 throughout the 
reasoning, so the demonstration is valid whether or not 
the positions of additional nodes are optimization param- 
eters. 

It must be mentioned that the absence of loops in 
the least dissipative network has already been specu- 
lated, but not formally proved, in the particular case 
of a constrained total pipe volume 0. Besides, in or- 
der to make connection with recent studies on networks 
minimizing the cost function X)fe w fcKfe| 7 0i 13 > ^ 
us rewrite the dissipation rate of a network with op- 
timized cross-sections as U — (kI) 5 J\ Ik |ifc| 7 ; where 

6 = 2m/ (m + n) and 7 = 2n/ (m + n) < 1 (this expres- 
sion was obtained by using Eqs. |2]and|2J). Banavar et 
al. 0] showed that the flow rate distribution satisfying 
the conservation law at every node, source and sink, and 
minimizing the cost function Wk |«fc| 7 is loopless when 

7 < 1. In that study however, the whole network struc- 



ture (i.e.: its topology and geometry) and the weight Wk 
of every link are set. The only optimization parameters 
are the flow rates ik (they do not necessarily obey to 
Ohm's law). Xue et al. |l2j obtained a similar result 
with the cost function ^2 k h |*fc| 7 , although the topol- 
ogy of the network is not necessarily settled, since the 
number and location of additional nodes are free to vary. 
However, this cost function differs from U, since k is not 
a constant, but depends on the geometry and topology 
of the network. Therefore, our result on the absence of 
loops cannot be deduced from those previous studies. 




(a) (b) 

FIG. 2: (a): the two adjoining pipes of a two-fold junction 
carry flows in opposite directions in order to satisfy flow rate 
conservation, (b): we can favorably replace the two adjoining 
pipes with a straight one: since the total channel length is 
shortened, the dissipation rate will be decreased for a fixed 
value of C n . 



Let us now show that the number of additional nodes 
is at most N — 2 in the optimal network, where N is the 
total number of sources plus sinks. This result limits the 
number of possible topologies for the optimal network. 
Note first that a loopless network (or tree) has one more 
node than it has links. So, the number of links in a net- 
work with A additional nodes is N+A—l. Since each link 
has two ends, the number of "incident lines" , summed 
over all the nodes, is 2 (N + A — 1). This number can 
be evaluated differently: let N p be the number of sources 
or sinks with p incident lines. Since each source or sink 
is linked to the rest of the tree, the smallest value of p 
for which N p has a nonzero value is 1, so J2 P >i N p = N. 
Similarly, let A p be the number of additional nodes with 
p incident lines. By definition, two-fold junctions can ex- 
ist only if its two links are not parallel. Such junctions 
cannot exist whenever their locations are optimization 
parameters (since condition 0] must be satisfied) . Two- 
fold junctions could a priori exist if their positions cannot 
be freely adjusted. However, they can be favorably (i.e.: 
with no increase of U and C n ) removed and their two 
adjoining pipes replaced with a straight one, the way it 
is depicted in Fig. [21 Thus, the smallest value of p for 
which A p is not zero is p = 3 in both cases, and the to- 
tal number of incident lines is: J2 P >iP^p + H P >3PAp- 
Comparing these two expressions for the number of in- 
cident lines, and considering that X)p>iP-^p — N and 
Zip>3 P^p — h appears that: 

A < N - 2, (9) 
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which was to be proven. 

Because of the broad definition of the flow resistance 
(Eq. [JJ , the results presented in this letter can be applied 
in various situations. For example, the m = 1 case can 
correspond to electrical current in wires, liquid flow in 
porous conducts, mass or heat diffusion in bars (provided 
that for the latter the bar lateral surface is insulated). 
The in = 2 case corresponds to the laminar Poiseuille 
flow in hollow pipes. Minimization can be done under a 
fixed lateral surface area (n — 1/2) if one wants to save 
the material required to build the hollow pipes, or a fixed 
volume (n = 1), if one wants to preserve the amount of 
liquid flowing through the network. 

Unfortunately, the results presented in this letter do 
not give insights into the way of building the optimal net- 
work practically, or even into the uniqueness of such an 
optimal network. In fact, as for the Steiner tree problem - 
which consists in finding the tree of minimal length inter- 
connecting a set of given points- this problem is likely to 
be NP-hard, meaning that the solution cannot be found 
without an exhaustive search of all the possible topolo- 
gies. However, the NP-hardness does not exclude the 
possibility of establishing basic properties on the geome- 
try and topology of Steiner trees [lj| • Similarly, we were 
able to address features on the structure of pipe networks 
minimizing the total dissipation rate under a global con- 
straint. Specifically, the upper bound on the number of 
additional nodes restricts the number of possible topolo- 
gies for the optimal network(s). These results make pos- 
sible the conception of efficient algorithms for computing 
the optimal pipe network problem [l2T |. In many situ- 
ations however, and especially for vascular systems of 
organisms, the capacity of the network to resist random 
injuries may also play a key role in its design. Obvi- 
ously, a reticulate network containing redundant paths is 
more adapted than an arborescent one for that purpose. 
Therefore, it is sometimes essential to look for a compro- 
mise between optimization of flow and robustness of the 
network. Definitely more works need to be conducted in 
that direction. 

B. Abou, S. Durand, and A. Rabodzey are deeply ac- 
knowledged for useful discussions and inspection of this 
work. 

APPENDIX 

In this appendix we prove that conditions [21 and 0] can- 
not lead to a local maximum of U . It is more convenient 
to show first that for a fixed value of U, a small devi- 
ation from condition 01 can lead to an increase of C n . 
Suppose that the location of a given node in the net- 
work is moved by a small amount Sr. The new length 
of each of the adjoining links can be easily evaluated 

/ 2 \l/2 

as: l' k = l k I 1 - 2e fe ■ Sr/l k + (5r/l k ) (where l k is 



its original length). Now let us adapt the cross-sectional 
areas such that each individual resistance is unchanged: 
r' k = r k . Thus, the dissipation rate is also unchanged, 
while the value of the total volume (or surface area) is 
now C' n — s »/^ OT+n " m y7™/ m j which can be expanded 
as: 

C' n = C n - 1 ^^2ste k .ST (10) 

j X ' 

Hence, when conditionals satisfied, we obtain that C' n > 
C n (with the assumption that m > n). We can then 
decrease C' n down to C n by reducing some cross-sectional 
areas. This can only imply an increase in U . Thus, we 
found a small perturbation from the extremum defined 
by Eqs. |21 and 01 leading to an increase in U for a fixed 
value of C n . We conclude that this extremum cannot be 
a local maximum of U . 
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